Origin of Multikinks in Dispersive Nonlinear Systems 
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We develop the first analytical theory of multikinks for strongly dispersive nonlinear systems, 
considering the examples of the weakly discrete sine-Gordon model and the generalized Frenkel- 
Kontorova model with a piecewise parabolic potential. We reveal that there are no 27r-kinks for 
this model, but there exist discrete sets of 27rA^-kinks for all iV > 1. We also show their bifurcation 
structure in driven damped systems. 



Non-equilibrium dynamics of many physical systems 
can be characterized by the creation and motion of topo- 
logical excitations or defects. In particular, when a 
nonlinear system possesses a degeneracy of its ground 
state, such excitations are kinks, the simplest and prob- 
ably most studied nonlinear modes. The concept of 
kinks is vital for many physical problems such as disloca- 
tion and mass transport in solids, charge-density waves, 
commensurable-incommensurable phase transitions, con- 
ductivity, tribology, Josephson transmission lines, etc [Q. 

In application to problems in solid state physics, the 
kink's motion is strongly affected by the inherent lat- 
tice discreteness. Earlier numerical simulations |^ of 
the kink's motion in a lattice described by the discrete 
sine-Gordon (SG) equation, also known as the Frenkel- 
Kontorova (FK) model demonstrated a number of 
interesting features not observed in the dynamics of soli- 
tons of integrable (both continuous and discrete) models. 
In particular, Peyrard and Kruskal found that a single 
kink becomes unstable when it moves in a discrete lat- 
tice at sufficiently large velocity, whereas two (or more) 
kinks are stable and propagate as multikinks. The former 
effect is associated with resonant interaction between a 
kink and radiation ||] , and resonances are even observed 
experimentally [Q. In contrast, the latter phenomenon, 
i.e. the formation of multikinks, "... has no clear ana- 
lytical explanation yet" (see p. 25). 

Recently, different physical systems have been studied 
numerically where multikinks are found to play an im- 
portant role. For example, multikinks are responsible for 
a mobility hysteresis in a damped driven commensurable 
chain of atoms ||^. In arrays of Josephson junctions, in- 
stabilities of fast kinks lead to the generation of hunched 
fluxon states also described by multikink modes j^. 

The main purpose of this paper is to provide the 
first step towards an analytical theory of multikinks in 
strongly dispersive nonlinear nonintegrable systems, in- 
cluding the analysis of the existence and codimension of 
A'^— kink states. In particular, we consider a weakly dis- 
crete SG model and demonstrate the existence of a finite 
number of multikinks due to a higher-order dispersion. 
We also find analytical solutions for multikinks and de- 



scribe the effect of an external field and damping on their 
existence and qualitative features. 

We consider the dynamics of a commensurable chain of 
atoms in a periodic substrate potential. In a normalized 
form, the equations of motion for the atomic displace- 
ments Un can be written as 

= V^(it(aO+'"n+l-W«)-V^nt(ao+Wn-'«ri-l)-W'sub("«), 

where Vint(w) is an effective interaction potential with 
the equilibrium distance oq, and Wsub(w) is a substrate 
potential with period a. For small anharmonicity, i.e. 
when \un+i — Wn| ^ do, the potential T4nt(w) can be ex- 
panded into a Taylor series to yield (see details in Ref. 
l|l|): iln- g{un+i + Un-i - 2u„) -|- W!.^^{un) = 0, whcrc 
g = V5J(^(ao). In the quasi-continuum limit, taking into 
account a higher-order dispersion, we obtain the normal- 
ized equation 

Utt - Uxx - Pu^xxx + W^s'ub(^*) = 0> (1) 

where W{u) has rescaled period 27r and, for harmonic 
interaction, (3 — a^/I2. 

Equation (|^) takes into account the effect of lattice 
discreteness through a fourth-order dispersion term, and 
for /3 = and W^^y^{u) = sin u, it transforms into the 
well-known exactly integrable SG equation that has an 
analytical solution for a single 2-K-kink moving with ve- 
locity V, u = 4tan~^{exp[(a; — vt)/^l — f^]}. Similar 
kinks exist for a rather general topology of the substrate 
potential Wsub(u) [Q. However, our aim in this paper is 
to study a new class of localised solutions of Eq. (|l|) for 
/3 ^ in the form of 27riV- multikinks for > 1. 

First, following the original study of Peyrard and 
Kruskal Q , we consider the harmonic substrate potential 

Wsub(w) = 1 - cos u. (2) 

We look for kink-type localised solutions of Eq. (|^) 
that move with velocity v (v^ < 1), i.e. we assume 
u{z) = u{x — vt). Linearizing Eq. (Q) and taking 
u{z) ^ e^^' , we find eigenvalues A of the form, 

= - 1) ± V(l 4/3], 
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so that for /3 > there always exist two real and two 
purely imaginary eigenvalues. Thus, the origin u = 
is a saddle-centre point and hence kinks, which are ho- 
moclinic solutions to u = (mod 27r), should occur for 
isolated values of v for fixed /3 (see Refs. 0,^). That is 
they are of codimension one. Moreover, this codimcnsion 
is only true if the solutions are themselves reversible, that 
is invariant under one of the transformations: 



i?i : u(mod 27r) ^ — u(mod27r), u" 



t 



-t, 



i?2 : u' -u\ u'" -m'", t ~> ~t, 
where prime stands for differentiation with respect to z. 
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FIG. 1. The four 47r-kink solutions of Eq. 
^sub(^) = sin u propagating at the given velocities 
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To find all solutions of this type, first we fix /? = 1/12, 
which corresponds to a = 1. Then, we perform numeri- 
cal shooting on the ordinary differential equation for u{z) 
using a well-established Newton-type method for homo- 
clinic/heteroclinic trajectories in reversible systems (see 
[^). The first result is that there exists no 2t: -kink so- 
lution at all, except in the artificial limit — > — oo (see 
comment below). Instead, we find a discrete family of 
47r- kinks; specifically there exist only four such solutions 
at four different values of v. The first solution has an 
analytical form |10[| 



u{z) 



8 tan ^ exp 



{(3/3)1/^} 



(3) 



where 



^47r 



= 1 - 2^/l3/i, i.e. for 
2/3. Other values are 



(2) 



our choice of /3, 
0.59498 . . ., 



0.42373 . . ., and v 



(4) 



0.21109. 



All these solu- 



tions are presented in Fig. 1. We may regard this discrete 
family as part of an infinite sequence of bound-states of 
two 27r-kinks that converges to the limit of infinite sepa- 
ration at a value of < 0. Actually the key parameter is 
fjL = 1 — , and further numerical evidence reveals that 
the bound states converge to /i — cx) at which value a 
27r-kink exists only formally. 

In addition to the 47r- kinks, numerics further reveals 
?;-values at which 27rA^-kinks occur for all iV > 2. Fig- 
ure 2 shows several examples of Gtt- and Stt- kinks. Ac- 
cording to a dynamical systems theory result H, on 



the existence of bound states of homoclinic solutions 
to saddle-centre equilibria in reversible Hamiltonian sys- 
tems, again thinking of the 47r-kinks as bound states of 
27r-kinks, one should expect to see precisely two Gir-kinks 
for each Air-kink. These would occur at fg^'^ satisfying 
v'^^l < v[^l < Wg'j^; all eight of which are depicted in Fig. 
2(a). Moreover, there would be two infinite sequences of 



SiT-kinks at v 



(ij)± 



871- 



such that z) 



from below as 



J 



oo and v. 



"4^ 



from above. Our numerical 



simulations have revealed precisely this structure of all 
multikink families. 
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FIG. 2. Examples of (a) &-k- and (b) Syr-kinks of Eq. 
with the potential (m) at given velocities. 




FIG. 3. Two-parameter continuation of 47r-kink solutions 
for the model (|l ) with the harmonic nonlinearity ^ . 



Finally, it appears that the above structure is largely 
independent of [3. Figure 3 shows the results of continu- 
ation (using the method for homo/heteroclinic orbits in 
the boundary- value software AUTO |jll|]) of the four 47r- 
kinks in the {v,\/\/]3) plane. These curves are almost 
identical to those obtained numerically in Ref. for the 
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discrete SG equation. Note that no curve passes through 
V = 1, they only reach there asymptotically as /3 — > 0. In 
the process the slope of each kink at its midpoint steep- 
ens, so that the solution becomes singular in the limit. 

It is important that the above numerical results may 
be verified by the construction of exact solutions in closed 
form when the substrate potential is approximated by a 
piecewise parabolic potential that generates in Eq. (|^) 
the effective force, W^^y^{u) = 

( u - 2mT : {2n - 1)tt + tt/2 < u < 2mT + n/2, 

\ {2n +1)7: ~u : 2mr + 7r/2 < u < (2n + l)7r + n/2. 

For simplicity we fix /3 = 1/12 and then look for 
kinks moving with velocity v — 1 — u^), by solving 
the piecewise-linear equation for u{z). This defines a 
four-dimensional dynamical system in the phase space 
{u,u',u",u"') e (-7r/2, 37r/2] x R^. The phase space is 
separated into two distinct domains: 

Region 1 : \u\ < 7r/2, Region 2 : 7r/2 < u < 37r/2. 

47r- kinks can be constructed by first noticing that, in or- 
der to be of codimension-one (i.e. occur at isolated v- 
values), they should be reversible under the transforma- 
tion Ri above. Since we can always translate by multiples 
of 27r, we look for solutions which satisfy, for some un- 
known Z2, the conditions: u(— 00) — > 0, 14(^2) = 27r, and 
'u"{z2) = 0, so that u{z) is in Region 1 for all z < 0, in 
Region 2 for < z < zi, for some unknown zi < Z2, and 
is in Region 1 again for all zi < z < Z2. 

The boundary condition can be satisfied by noticing 
that such solutions at z = (the first point of tran- 
sition between Regions 1 and 2) satisfy u{0) — tt/2, 
w'(0) = A7r/2, u"{0) ^ A^7r/2, and u"'(0) = X^Tr/2, 
where = 6{^yJt? + 2/3^^ — ^) is the unique real positive 
eigenvalue of the linear system in Region 1. Hence the 
asymptotic boundary condition ai z — —00 in Region 1 
becomes an initial condition at z = for u in Region 2. 
The general solutions in Regions 1 and 2 are: 

ui{z) = Aie^'^ + Bie^^ + Ci cos(wz) + Di sin(wz) 

and, providing /i > /imin :— ^2/37r, 

lt2(2:) = A2 C0s(tJiZ + B2) +C2 C0S(UJ2Z + D2) 

where — 6(y7?+2/37r + /i) and a>^_2 = 6(^ ± 
^A*^ - 2/37r), Aj, Bj, Cj, and Dj are unknown coef- 
ficients. Therefore, we can explicitly solve for the co- 
efficients to find U2{z) in closed form. This expression 
defines an implicit equation for 21; ^2(2:1) — 37r/2. The 
value of ^2(2:1) and its derivatives then defines initial con- 
ditions at z = zi, hence determining the constants Ai, 
-Bi, Ci, and Di. This in turn defines Z2 implicitly as 
Wi(z2) = 2'K. To have a 47r-kink we additionally require 
w'l (Z2) = 0, and so should only expect to find zeros of this 



final quantity by varying v. Hence we can define a 'test 
function' for 47r-kinks K{v;zi,Z2) ■= U2{z2)- Using the 
above construction, this K can be written in closed form 
in terms of w, zi, and Z2. The unknown transition points 
zi_2 are the solution to given transcendental equations, 
in each case only the first solution of which has meaning. 




FIG. 4. The function K{v) for multikinks of Eq. (0) 
with piecewise parabolic potential. Inset: the corresponding 
47r-kinks and their associated u- values. 

Figure 4 shows a graph of if as a function of /i = 1 — 

G (/^min: l)i which has been computed using MAPLE 
with the implicit equations solved for their smallest pos- 
itive solutions. The five zeros of K correspond to 47r- 
kinks, graphs of which arc shown in the insert to the fig- 
ure. These zeros occur for v = 0.64064609, 0.49870155, 
0.37835717, 0.26634472, 0.14477294. It is also possible 
to construct solutions for /i < /imin in a analogous man- 
ner, but with the solution in Region 2 replaced by one 
corresponding to complex eigenvalues. This gives the ad- 
ditional solution for v — 0.833706. 

In this way, we find analytically a finite set of v- 
values giving 47r-kinks for the piecewise parabolic poten- 
tial model, having qualitatively the same structure as the 
solutions found numerically for the sinusoidal nonlincar- 
ity (||). One could go on to construct 27rA^-kinks for 
N > 2, but the calculations presented already serve to 
corroborate the earlier numerical results. 

To complete the analysis of the kinks, we would like 
to mention that the short-wavelength instability of the 
nonstationary continuous model (|l|) due to the term 
be easily removed by introducing an equiva- 
lent higher-order dispersion via a mixed derivative term 
00. 

To analyse the robustness of multikinks in realistic 
physical systems, we add to the right-hand side of Eq. 
(0) the driven damped term F — 6ut, where F is an ex- 
ternal DC force and i5 is a damping coefficient (see, e.g.. 
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g). Importantly, for each of the kinks so far found, it 
is possible to use numerical continuation to trace curves 
that lie on sheets in the parameter space {v,S,F) cor- 
responding to the existence of multikinks. For example, 
taking the explicit 47r-kink solution given by Eq. (H), a 
curve was computed a.t f3 — 1/12 in the {F, v(5)-plane with 
fixed \v\ — y^2/3, reaching a maximum with respect to 
S at \v\S = 0.069326. Taking the fixed value \v\S = 0.05 
from this curve the locus of kinks in the (w, -F)-plane can 
then be traced out, as depicted in Fig. 5. 
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FIG. 5. The kink's velocity v < and maximum ampli- 
tude against external force F > for the simplest 4n kink 
with v5 — 0.05. Note that each limb of the curve spirals back 
on itself. The insets illustrate example solutions on the locus. 

Three interesting features can be noted from this curve. 
First, all kinks have developed oscillations around the 
equilibrium close to u = 47r. This is because, for 5 > 0, 
the corresponding equation for travelling waves is no 
longer Hamiltonian or reversible, and the linearization 
around the asymptotic value — sin^^ F now has three 
stable eigenvalues, two of which have non-zero imaginary 
part. These oscillations may be regarded as radiation 
that travels at the kink 's velocity, as was earlier observed 
in direct numerical simulations [|6|. Second, v and F have 
opposite sign for these results. When F and v have the 
same sign, only kinks with non- decaying oscillations in 
the tails can be found. Third, note that the computed 
curve ends at a point where a transition takes place in- 
volving a heteroclinic connection with u w Stt. This sug- 
gests that TT-kinks are possible for sufficiently large F. 

Finally, we mention that the case /? < in Eq. (|l|) 
can also occur in generalised nonlinear lattices provided 
we take into account the next-neighbor interactions, e.g. 
due to the so-called helicoidal terms in nonlinear mod- 
els of DNA dynamics In this case, the analysis is 
much simpler and, similar to the nonlocal SG equations 



p4[ , leads to the continuous families of multikinks pa- 
rameterized by V. From the mathematical point of view, 
for (3 < the origin changes from a saddle-center to a 
saddle-focus, and rigorous variational principles | [T5[ | give 
families of stable 27rA^-kinks for all iV > 1. 

In conclusion, we have developed the first analyti- 
cal theory of multikinks in strongly dispersive nonlin- 
ear systems, considering the important examples of the 
generalized FK model with the sinusoidal and piecewise 
parabolic potentials. We have revealed, numerically and 
analytically, the existence of discrete sets of 27rA^- kinks. 
We believe that general features of multikinks and the 
physical mechanism for their formation are similar in 
many other strongly dispersive nonlinear models. 
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